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Abstract 

•l-H 

43 We present robust radio frequency (rf) pulses that tolerate a factor of six inho- 

, Ph 

mogeneity in the Bi field, significantly enhancing the potential of toroid cavity 

J> resonators for NMR spectroscopic applications. Both point-to-point (PP) and uni- 

00 

tary rotation (UR) pulses were optimized for excitation, inversion, and refocusing 

(N 

using the gradient ascent pulse engineering (GRAPE) algorithm based on optimal 
control theory. In addition, the optimized parameterization (OP) algorithm applied 

o 

t-h to the adiabatic BIR4 UR pulse scheme enabled ultra-short (50/iis) pulses with ac- 

• i-H ceptable performance compared to standard implementations. OP also discovered a 

new class of non-adiabatic pulse shapes with improved performance within the BIR4 
framework. However, none of the 0P-BIR4 pulses are competitive with the more 
generally optimized UR pulses. The advantages of the new pulses are demonstrated 
in simulations and experiments. In particular, the DQF COSY result presented here 
represents the first implementation of 2D NMR spectroscopy using a toroid probe. 
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1 Introduction 



Practical NMR applications require pulses that provide robust performance with respect 
to experimental limitations, such as resonance offset effects and rf inhomogeneity. A 
simple rectangular pulse delivered with a perfectly homogeneous rf-amplitude of 20 kHz 
(pulse length 12.5 /is) transforms M z to M x with a fidelity of 99% over an offset range of 
only ±2.8 kHz. Miscalibration or inhomogeneity of the Bi field exceeding ±9% reduces 
M x on resonance below the desired fidelity. 

Increased tolerance to offset and/or rf inhomogeneity can be achieved using composite 
(JU |2]) and shaped pulses. More recently, efficient pulse design using optimal control 
theory (BJ has made it possible to establish physical limits to pulse performance (0 El [7J EJ 
[91 HOI HH [121 H3l [HD . Short (100 /is), broadband, PP excitation pulses (98% fidelity) have 
been optimized that achieve bandwidth to peak rf ratios of 2 with a ± 10% tolerance to rf 
inhomogeneity applicable to modern high resolution NMR probes dE}. In an extreme case 
study of dual compensation, a 1 ms pulse was found that provides excellent broadband 
excitation over an offset range of 50 kHz for miscalibration of the Bx field anywhere in the 
range 10-20 kHz (|10j) . Much larger spatial variations of the Bi field limit appplications 
involving, e.g., surface coils, ex situ NMR, and high field imaging. 

Here we focus on toroid cavity probes (1T5| IT6| H7I) . which can be built to tolerate high rf 
amplitude, high pressure, and high temperature, but at the price of large rf field inhomo- 
geneity. A toroid cavity detector is a unique NMR resonator with a wide range of potential 
applications for in situ reaction studies at high-pressure and/or high-temperature (jTBl . 
Whereas the defined B\ field gradient of the coil is an advantage for some applications 
(for example, rotating frame NMR imaging on the micrometer scale (19)), it is a serious 
problem in spectroscopic applications, severely limiting the potential of toroid NMR. For 
the toroid probe designs considered here, the ratio between the minimum and maximum 
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Bi field in the sample volume is a factor of six. Other detector geometries produce even 
larger inhomogeneities ^TTJ). 

We present efficient point-to-point (PP) and unitary rotation (UR) pulses optimized for 
toroid probe applications using the GRAPE algorithm (20). PP pulses rotate one specified 
initial state about a fixed axis to a desired target state, whereas UR pulses transform any 
orientation of the initial magnetization about the same axis. While composite PP pulses 
developed specifically for toroid detectors exist, (I2"2"j) . no UR pulses, which are crucial for 
multi-dimensional NMR, have been developed previously for these probes. The perfor- 
mance of individual toroid pulses is characterized both theoretically and experimentally 
in section 3. We demonstrate the performance and relevance of the new pulses for toroid 
probe NMR in a DQF COSY experiment, where significant gains in signal amplitude are 
found compared to experiments using conventional pulses. 



2 Toroid probes 



The first use of toroid coils in NMR spectroscopy was reported in 1983 ([TBI 1231) • Toroid 
rf coils and toroid cavities provide excellent signal-to-noise ratio compared to conven- 
tional Helmholtz or saddle coils, and very high rf amplitudes on the order of 100 kHz 
can be reached. In toroid cavity autoclaves (1T8"| |2~4"|) . pressures of up to 300 bar and 
temperatures of up to 250 °C are attainable, making it possible, e.g., to investigate in 
situ reaction dynamics of homogeneously catalyzed reactions, such as the cobalt-catalyzed 
hydroformylation process (|2"5"1 |2T| fI7\ |2"8"]) . 

In a typical toroid probe with cylindrical symmetry, the sample is located between a 
minimum radius r OT j n (given by the radius of the central conductor) and a maximum radius 
fmax- The B\ field for a toroidal geometry varies inversely with radial distance r as detailed 
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more fully in (ITT j l29 | l30l). Hence, the corresponding rf amplitude, v r f(r) = r yB 1 (r)/(27r), 
can be expressed relative to its smallest value at r max as v r f{ r ) ^ v r f{r max ) jr to obtain 



The time dependence of the rf pulse can then be written simply in terms of an amplitude 
modulation function < ait) < 1 as 



applied to the rf spatial profile at any position r in the toroid. 

For the toroid probe used here, r m j n = 1 mm and r max = 6 mm. A conventional rectan- 
gular pulse thus rotates spins near the central conductor by a flip angle that is six times 
larger than the flip angle experienced by spins near the outer wall of the toroid resonator. 
This large and well defined rf inhomogeneity of toroid probes can be exploited, e.g., in 
spatially resolved diffusion measurements and imaging (TIT H IHT | 133]) . However, the large 
rf inhomogeneity has limited spectroscopic applications of toroid probes to relatively sim- 
ple experiments. In the following, we remove this limitation by developing pulses with 
the necessary high tolerance to rf inhomogeneity. 

3 Pulse optimizations and applications 

The GRAPE algorithm for pulse optimization is discussed in detail in the cited references 
on optimal control. A quality factor, <3>, for pulse performance is defined which, in turn, 
provides an efficiently calculated gradient for iterative improvement of pulse performance. 
Most generally, the quality factor is a quantitative comparison between the state of the 
system and some desired target state. The gradient therefore also depends on the system 
state — in the present case, the magnetization M. Modifications in the basic algorithm 




(1) 




(2) 
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that are required for the large and well-defined rf spatial inhomogeneity in toroid probes 
requires some elaboration. 

3.1 Simulating pulse performance in a toroid 

For a general rf pulse, each combination of offset, v ff, and rf amplitude, v r f(r), produces 
a potentially different transformation of the initial magnetization. The goal of pulse 
optimization is to find a particular rf pulse that produces the same transformation for all 
the scalings of the rf amplitude due to spatial inhomogeneity and all the desired offsets. 
A gradient giving the proportional adjustment to make in the rf pulse components to 
improve performance at each u Q ff and v r f{r) can be efficiently calculated for point-to- 
point (PP) pulses (0 [20]) an d f° r unitary rotation (UR) pulses ([2"U]) . The total gradient 
is obtained by averaging these constituent gradients over offset and rf inhomogeneity. 

For the small volumes and relatively small deviations from homogeneity seen in standard 
NMR probes, giving equal weight to different possible spatial values of the rf is sufficient 
to provide accurate simulations of pulse performance relative to experiment. In a toroid 
probe, the effect of the large and well-defined spatial inhomogeneity on both the spa- 
tial dependence of the transformed magnetization and the spatially dependent detection 
sensitivity must be considered. 

The signal from the toroid cavity depends on the total contribution from spins at each 
radius. Signal is proportional to the detection sensitivity per spin times the number of 
spins within a cylindrical sample slice of height h, inner radius r, and outer radius r + Sr. 
The slice volume 8V = 2irrh5r (and, thus, the corresponding number of spins in the slice) 
increases linearly with r. By the principle of reciprocity (l2Til . the detection sensitivity is 
proportional to v r f and hence proportional to 1/r (c.f. Eqs. 1 and 2). Therefore, the 
signal from a cylindrical slice volume 5V is independent of r, and we can define an effective 
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magnetization vector representing the sample in a toroid probe as 

M eff (t) = / M(r,t)dr (3) 

^max Train Jr m i n 

with equal weighting at each radius. The (x, y) components of the detected signal are 
proportional to the (x,y) components of M e -^(t), respectively. The magnetization vec- 
tor M(r, t) in Eq.([3]) resulting from the applied rf is calculated starting from thermal 
equilibrium M = (0, 0, 1) using the rf amplitude v r f(r) given in Eq.Q. 

In the numerical simulations, the integration in Eq.([3]) is approximated by a discrete 
sum. For each offset, the gradient (which depends on M) is averaged over the range of 
rf spatial variation by one of the two methods outlined in Appendix A. The first method 
samples M(r, t) at equally spaced r. Since v r f varies as 1/r, this has the effect of coarsely 
digitizing the rf for small values of r and sets the stepsize Ar required for accurate 
simulations of rf inhomogeneity in the toroid. At large r, however, Ar is more accurate 
than necessary. Alternatively, the second method samples M[^v/(r), t] at equally spaced 
u r f, multiplied by the proper weight for each rf frequency that accurately represents the 
nonlinearity of u r f(r), as derived in Appendix A. The gradients resulting from either 
procedure, averaged over the range of rf spatial variation for each individual offset, are 
subsequently averaged with equal weight over the range of offsets u a fj to give the overall 
gradient for the performance factor $. 

Optimizations for a range of experimental parameters obtained efficient and robust pulses 
for ratios r max /r m i n as high as 100, which is significant for extending the present results to 
additional toroid applications. An experimental upper limit v™j ax = 25 kHz on v r j[r max ) 
for our toroid cavity resonator (r min = 1 mm, r max = 6 mm) limited the time-dependent 
pulse amplitude at r max to the range < v r f{r max , t) < v™ x in the optimizations. We 
consider offsets v g in the range ±1.5 kHz, corresponding to a 15 ppm : H chemical shift at 
a spectrometer frequency of 200 MHz. Different experimental settings are accommodated 



7 



by expressing frequencies relative to the limit u^j ax . In these relative frequency units, the 
current experimental setting corresponds to an offset range of —0.06 < f ff/ 1/ ^f' x < 0.06. 
Pulses were optimized starting from a set of random initial pulses. Each pulse was digitized 
in time steps of 0.5 fis duration and, at each time step tj, the rf amplitude modulation 
function a(tj) and the phase 4>(tj) were optimized. 



3.2 Point-to-point (PP) pulses 

Point-to-point transformations rotate one specified initial state about a fixed axis to a 
desired target state. For example, a 90° PP pulse for the rotation M z — > M x will not 
in general rotate any other component 90° about the ?/-axis. An optimal hard pulse 
achieves a maximum fidelity for this PP transformation of only 74% (jTTJ) (see Appendix 
B). A conservative lower bound for the required duration of a composite or shaped PP 
excitation pulse to achieve excellent overall performance can be estimated in comparison 
to the duration of a 90° hard pulse with amplitude set to the toroid upper limit i>^j ax = 
25 kHz available at r max . This results in a conservative estimate for the lower bound of 
l/(4z/j a:c ) = 10 fis. The minimum time to reach a fidelity of, e.g., 99% is expected to be 
significantly longer than 10 /is, since this hard pulse would generate a 540° rotation at 
r min . Hence, additional time is needed to create uniform excitation for all rf amplitudes 
in the sample and all offsets of interest. 

To characterize the performance limits of robust excitation pulses, we considered pulse 
lengths T in the range 2.5-40 //s. Multiple optimizations were performed starting with 
random pulse shapes at each T to obtain the best overall efficiency for transforming 
magnetization M z — > M x in the toroid probe. The resulting time-optimal pulse (TOP) 
curve is shown in Fig. [I] The TOP curve reaches a value of 99% for the quality factor 

$ PP = M e J f (4) 
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Figure 1: Time-optimal pulse (TOP) curves showing optimal pulse performance as a function of pulse 
length for the PP transformation M z — > M x applied to the toroid geometry and resonance offset range 
described in the text. The optimal length on resonance hard pulse derived in Appendix B (T — 3.74 /xs) 
is plotted (open diamond) for comparison. M%ff generated by three optimized composite pulses based 
on the phase-modulated excitation pulses in (22) is shown as gray squares. An enlarged view of the 
region where Mfj.ff > 0.98 is shown in B. Optimized composite pulse parameters are provided in the 
Supplementary Material. 
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at T rs 22 /xs, which is only 2.2 times longer than the conservative estimate of 10 [is for 
the lower bound given above. 

Existing composite pulse schemes (1221) that are also robust to a broad range of rf inhomo- 
geneity were optimized for use in the toroid and are plotted (open squares) for comparison. 
Pulse parameters are provided in the Supplementary Material. For this experimental set- 
ting, a composite pulse constructed with an optimal duration of 25.62 /js achieves a quality 
factor of 99.30%. This pulse approaches the performance shown in the TOP curve of Fig. 
[TJ indicating that it is close to the performance limit for this pulse duration. However, a 
shorter optimized pulse of duration 22 /is is able to achieve the same quality factor as the 
composite pulse. 

More detailed performance of this constant amplitude, phase-modulated composite pulse 
is shown in Fig. [2] compared to an optimized phase-modulated pulse of the same length. 
There are many different optimized pulses that achieve performance similar to the pulse 
shown. This partcular pulse was chosen because the profile of its phase modulation is 
reminiscent of the composite pulse. Pulses with even more striking similarity to the 
composite pulse shape can be found at T = 29 /is with slightly improved performance. 

The close match that is possible to achieve between theoretical and experimental pulse 
performance is illustrated in Fig. [3] for a high-fidelity excitation pulse (T = 100 /is) that 
transforms M z — > M x with a quality factor of 0.9993. The value of M x is plotted as a 
function of resonance offset and radial variation of the rf field strength, scaled as described 
in Fig.[2j The region for which the pulse was optimized (—0.06 < ^o///^ ax — 0.06 and 
1 — v rf{f) I v™f ax < 6 ) is indicated by a box in the figure. The experiments were performed 
on a Bruker AV 250 spectrometer equipped with a conventional 5mm QNP probe. A 
Shigemi tube was filled to a height of 4 mm with 1% H 2 in D 2 doped with copper 
sulfate to enhance relaxation and reduce experimental time. In this experimental setting, 
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Figure 2: Comparison of the composite excitation pulse from Fig. [T] with T = 25.62 /zs designed 
according to Ref. (|2"2")) and the optimzed pulse of equal length. Both pulses are purely phase-modulated 
and transform M z — > M x . (A) Composite pulse phase. (B) Excited magnetization plotted as a function 
of resonance offset v $ f and radial variation of the rf field strength v r f (r) given by Eq. ([!]) . Both axes 
are normalized to the maximum rf field v^ ax — 25 kHz available at r max in the toroid probe. (C) and 
(D) Corresponding plots for the optimized pulse. The composite pulse performance is very close to the 
optimized pulse, but, as shown in Fig.[T] a shorter optimized pulse of length 22 fis can be used to achieve 
performance similar to the composite pulse. 
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Figure 3: Simulated and experimental performance of a high-quality PP pulse (T = 100 (is) trans- 
forming M z — > M x plotted as a function of resonance offset and radial variation of the rf field strength, 
as described in Fig.[2j The value of M x in the region for which the pulses were optimized (indicated by 
a black rectangle) is 0.999. 

the rf inhomogeneity is on the order of 1% and can be neglected. We systematically 
scanned offset (50 steps) and the rf (50 steps) in a series of 2500 experiments. An excellent 
match was found between the simulated and experimental performance of the pulse. 

The TOP curve thus provides a useful benchmark for past and future PP excitation pulses 
to be used in toroid probes. Although improved PP excitation pulses were found in the 
current setting for 1 H applications of toroid probes, the improvement in excitation perfor- 
mance compared to previously developed composite pulses was relatively small. However, 
this improvement provides the option for decreasing pulse lengths by 10-15%. More im- 
portantly, the methods presented here can be applied to applications where larger offset 
ranges need to be covered, e.g., for 13 C spectroscopy or higher spectrometer frequencies, 
where even greater performance gains compared to existing PP pulses can be expected. 

3.3 Universal rotation (UR) pulses 

Universal rotation pulses that transform any orientation of the initial magnetization about 
the same axis are desirable whenever more than one magnetization component has to be 
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rotated in a controlled way. This is the case for most multi-dimensional NMR experiments, 
which cannot be efficiently implemented based on PP pulses alone. In the following, we 
characterize the performance of UR pulses suitable for use in the toroid probe, analogous 
to our treatment of PP pulses. The primary difference is that our quality factor for 
assessing performance becomes (1201 

$ UR = Re(U target | U eff (T)), (5) 

where U target is the desired (ideal) propagator for the transformation, U e ff(T) is the actual 
effective propagator at the end of the pulse, and the trace operator Tr sums the diagonal 
elements of the resulting matrix product. This is equivalent to combining (averaging) the 
quality factors for separate PP transformations about the coordinate axes. For example, 
a 90° universal rotation of single-spin magnetization about the y-axis axis averages $pp 
for M z -»■ M x , M x -M z , and M y ->> M y . 

Figure [4] shows the TOP curves we found for the performance limits of UR 90° and UR 
180° pulses suitable for use in the toroid. These pulses require approximately twice the 
time of PP pulses to approximate an ideal quality factor of one. This is consistent with 
earlier work on the construction of UR pulses from two PP pulses (1341) . The composite 
pulses of Ref. (l22l) are not able to achieve the UR performance of a simple hard pulse in 
the toroid. 

To date, the only other UR pulses with sufficient tolerance to rf scaling and resonance 
offset that might make them suitable for use in toroid probes are the adiabatic BIR-4 
pulses (j35j) . However, pulse lengths of ~ 50 fis for the UR pulses obtained here, requiring 
an extremely rapid 12.5 /as half-passage in each BIR-4 segment, are too short for good 
adiabatic performance. The peak performance we were able to obtain from BIR-4 at one 
particular combination of offset and rf scale factor was $>ur — 0.75. The performance over 
the volume of the toroid is considerably less. The BIR-4 pulse was constructed using the 
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Figure 4: Time-optimal pulse (TOP) curves showing optimal performance as a function of pulse length 
for the universal rotations UR 90° (A) and UR 180° (B) applied to the toroid geometry and resonance 
offset range described in the text. By comparison, the quality factor for the best optimally parameterized 
([36)) OP-BIR4 pulse we were able to generate for T — 50 [is is 0.9588. The performance of three optimized 
UR 180° pulses constructed by the procedure of Ref.(fTT[) from phase- modulated excitation pulses in 
is indicated by circles in (B). Optimized composite pulse parameters are provided in the Supplementary 
Material. 

tanh/tan functions for amplitude/frequency modulation ( |35|) . For a given performance 
level, they allow shorter pulse lengths than other adiabatic modulation schemes (1371) . 
Standard values of tan/t = 20 and £ = 10 or £ = 20 were used for the tanh/tan shape 
parameters. The frequency sweep was then adjusted for optimal performance. 

Thus, standard BIR-4 is unsuitable for such short pulse lengths. However, to more fully 
understand what is possible at T = 50 /is, u^j ax = 25 kHz, and the desired range of 
resonance offset and rf tolerance, we optimized the shape parameters and the frequency 
sweep for the tanh/tan pulse using the OP algorithm ( 136]) to obtain tanK = 17, £ = 71, 
and a frequency sweep of 602.5 kHz for each half-passage. The resulting quality factor, 
<&ur = 0.959 over the full volume of the toroid volume, is at the lower limit of the 
acceptable performance range considered in Fig.|4} However, this optimally parameterized 
OP-BIR4 pulse is not adiabatic and only retains the basic shape and form of BIR-4. More 
dramatic deviations from adiabaticity were found for the parameter values k = 10.35, 
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£ = 20.5, and frequency sweep 202.5 kHz, which produce the same &ur = 0.958. This 
large value of k, slightly greater than 6(7r/2), produces 1.5 "copies" of the standard 
BIR-4 frequency modulation scheme {n max = tt/2) in each half-passage segment, while 
the tanh amplitude modulation still has one cycle varying between and 1 in each segment. 
The limits of non-adiabatic performance using standard amplitude/frequency modulation 
functions is beyond the scope of the present paper. We simply note that a pulse length 
of 200 /is is needed before OP-BIR4 can achieve the same quality factor ($ur = 0.999) 
as the 50 /xs pulses derived for the TOP curve. 

One particular 50 /is 90° UR pulse from the TOP curve set is shown in Fig. [5] along with 
details of its performance throughout the optimization region. This pulse was chosen for il- 
lustration because it has interesting qualitative similarities to the BIR-4 amplitude/phase 
modulation. However, the frequency modulation profile of the optimized pulse (obtained 
as the derivative of the phase modulation, but not shown) makes it clear how different it 
is from a standard adiabatic frequency sweep. Figure|5] serves to illustrate the qualitative 
similarities that BIR-4 may have to an optimized pulse rather than the other way around. 
Adiabatic pulses are an elegant solution to the problem of large rf inhomogeneity, provid- 
ing a simple physical picture that is easy to understand. But optimal control shows there 
are many other solutions that are not so intuitive. 

3.4 Toroid NMR: Implementation of DQF COSY 

To demonstrate the performance of the optimized pulses in two-dimensional experiments, 
we chose the DQF COSY experiment (38) as a particularly simple illustrative example. 
Only the second and third pulses of this three-pulse sequence need to be UR pulses, 
but the same UR 90° pulse was used for all three pulses as a more stringent test of its 
capabilities. 
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Figure 5: The amplitude a(t) and phase tp(t) modulations of an optimally parameterized (36) OP- 
BIR4 pulse (see text) are shown in (A) compared to one particular 90° UR pulse of length 50 /us in (D) 
that achieves the performance limit represented in the TOP curve of Fig. [4] Neither optimized pulse is 
adiabatic, but the comparison illustrates the qualitative similarity the BIR-4 modulation scheme might 
have to a more generally optimized pulse. The effective rotation axis is plotted for various offset and 
rf scale factors below each pulse, (B) and (E), and the corresponding quality factor §ur = Re(f/go° | 
U e ff(T)} is plotted as a function of offset and rf scale factor in (C) and (F). The effective rotation axis 
for an ideal 90° UR pulse is the y-axis for all offsets and rf scales. For each (normalized) resonance offset 
v ff/v™f ax equal to —0.06 (red), (black) and 0.06 (blue), 21 effective axes with rf scale factor r max /r 
ranging from I to 6 were calculated. The effective rotation axes for pulses A and D show their deviation 
from the ideal, with further detail provided in the contour plots showing the quality factor achieved. In 
contrast to these two optimized pulses, BIR-4 constructed using standard parameter values (see text) 
gives <&ur < 0-75 over the desired range of offsets and rf scale factors. 
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Two sets of experiments were performed using a conventional probe for the first set and a 
toroid probe of dimensions r min = 1 mm and r max = 6 for the second. In each experiment, 
spectra were acquired using either conventional hard pulses for all pulses in the sequence 
or optimized pulses. 

The first set of DQF COSY experiments were recorded on a Bruker AV 600 spectrometer at 
600 MHz proton resonance frequency equipped with a conventional probe with negligible 
rf inhomogeneity. As a simple model system we used cytosine as a two spin system. The 
sample was prepared using a saturated solution of cytosine in a mixture of DMSO-de/D20 
(10/1, vol/vol) which was doped (139|) with paramagnetic Chromium(III)acetylacetonate 
to reduce the longitudinal relaxation time from 3.7 seconds, initially, to 0.15 seconds to 
accelerate experimental data acquisition. 

The upper limit for the rf amplitude in the conventional probe is 30 kHz, which scales to 
5 kHz to match the weakest rf amplitude in the toroid probe at r max = 6 mm. Since this 
is one fifth the value the pulse was originally optimized for, T must be increased by the 
same factor from 50 /is, originally, to 250 /is. The carrier frequency was offset 300 Hz from 
the N-Ci/=CH diagonal peak, corresponding to the maximum (positive) optimized offset 
range ^off/^f ax = ±0.06. The rf amplitude of the 90° hard pulse (u r f = 17.5kHz, T = 
14.29 fis) was chosen at the middle of the range of rf variation (5-30 kHz), corresponding 
to an rf scale factor of 3.5 times the minimum amplitude. 

Experimental results for the N-Ci/=CH diagonal peak are plotted as contours in Fig.[6]for 
rf scale factors r max /r equal to 1, 2.25, 3.5, 4.75 and 6. Simulated results are shown on the 
left. There is an excellent match between the simulated and experimental performance 
of the pulse. Although hard pulses work only within a very limited range of r max /r, the 
optimized UR pulse provides excellent performance over the entire rf spatial variation 
appropriate for the toroid. 
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Figure 6: Simulated (left) and experimental (right) DQF COSY spectra of cytosine. The diagonal peak 
of N-C-ff=CH is shown. Two series of spectra using rectangular and shaped pulses are displayed for 
rmax/r = 1.00, 2.25, 3.50, 4.75, 6.00 (c.f. text) 
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Figure 7: DQF COSY spectra of ethyl crotonate recorded on a toroid probe using (A) hard and (B) 
optimized UR pulses. Slices through the cross peak circled in each spectrum are displayed in (C) and 
(D), showing a signal gain of 2.2 for the experiment using the optimized UR pulse. Other peaks showed 
signal gains ranging from 1.7 to 2.7 relative to the hard pulse experiment. 
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The second set of experiments was performed on a 200 MHz Bruker Avance DRX spec- 
trometer using a homebuilt toroid probe with r min = 1 mm and r max = 6 mm. Details of 
the toroid design can be found in Ref. (l3ip . Use of an earlier generation console produced 
systematic errors in the pulse waveform if the pulse generator did not have sufficient time 
to respond to changes in pulse phase. To provide sufficient time, the pulse increment was 
first increased from 0.5 fis to 0.85 [is. Each pulse increment was preceded and followed 
by zero-amplitude sub-pulses of the same duration and phase to provide the necessary 
response time. The net affect was to increase T from 50/is to 357 /is. This scales the 
optimized peak rf and resonance offset range of the original 50/xs UR 90° pulse by the 
ratio of the pulse lengths. Instead of adjusting experimental settings to fit the modified 
pulse, we optimized a new 90° UR pulse for the same rf and resonance offset parameters 
considered throughout the present paper, but with a duration of T = 357 fis, a pulse 
increment duration of 0.85 fis, and alternating zero-amplitude sub-pulses as described 
above. 

Contour plots of complete 2D DFQ COSY spectra of ethyl crotonate obtained using a 
toroid probe are shown in Fig. [7] for rectangular (A) and optimized UR (B) pulses. The 
carrier frequency was set to ppm for both spectra, so that the chemical shift range for 
which the pulse was optimized ranges from —7.5 ppm < 5 < +7.5 ppm. Nevertheless, 
as the experimental results demonstrate, our UR pulse provides very good performance 
beyond the optimization range, up to +10 ppm. Simulations indicate an extended range 
of ±10.5 ppm in which the pulse provides acceptable performance. Slices from the spectra 
are displayed in C and D, highlighting the cross peak marked by circles in Fig. [7] A and B. 
For this cross peak, a signal gain of 2.2 is observed relative to the hard pulse experiment. 
For the other peaks of the spectra, 1.7 to 2.7-fold signal gains were obtained. Even larger 
gains can be expected for experiments with more pulses, such as the HMBC experiment. 
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4 Conclusion 



We have demonstrated the power of optimal control based methods for deriving rf pulses 
that tolerate the large B\ spatial inhomogeneity present in toroid probes. We considered 
both point-to-point (PP) rotations that transform a particular initial state to a desired 
final state and universal rotations (UR) that transform any initial state about a well- 
defined fixed axis to a desired final state. The limits of PP and UR pulse performance as 
a function of pulse length were provided for a fairly typical toroid probe geometry that 
produces a factor of six variation in the Bi field strength over the radial dimension of the 
toroid. Such pulses are crucial for realizing the full potential of toroid detectors for unique 
applications such as in situ reaction studies at high-pressure and/or high temperature. 

Optimized 90° UR pulses were utilized in a DQF COSY experiment to implement the 
first two-dimensional spectroscopic application using a toroid probe. Factors 1.7-2.7 gain 
in signal intensity were obtained compared to the conventional experiment using a hard 
pulse. Even larger gains can be expected for other homonuclear and heteronuclear exper- 
iments that employ more pulses or require performance over a larger range of resonance 
offsets. These kinds of experiments are important for future development of the field, but 
the present work is a significant step in opening the door to a wealth of new applications 
for toroid probes. 

In addition to PP excitation pulses, PP inversion pulses were also optimized and an ex- 
ample is given in the supplementary material. Excitation and inversion pulses considered 
in the text are available in electronic form at 

http: / / www.org. chemie.tu-muenchen.de/glaser/Downloads. html 

and can also be found in the supplementary material. 
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Appendix A: Discretizing methods for the calculation 
of the effective magnetization vector M e ^ 



The effective magnetization vector M e -^(t) is defined in Eq. (3) simply as an integral of 
the magnetization vector M(r, t) over the radius r from r m j n to r max . Here we consider 
two different methods to approximate the integral by a finite sum. 



Discretization method A based on equidistant steps Ar : 

The integration for M e -^(t) in Eq.([3]) can be approximated by sampling the radius r 
at N equidistant values r k {k = 1, 2, . . . , N) separated by Ar = {r max — r min )/N, with 
n = r min + Ar/2 and r N = r max - Ar/2, to give 

1 - 

k=l 

This uniform sampling in r, however, corresponds to a non-uniform sampling of the rf 
amplitude v r f{ r k) given by Eq.Q. For r min = 1 mm, r max = 6 mm, and z/j ax = 25 kHz, 
the rf frequency difference between u r f(l mm) = 150 kHz and v r f(1.5 mm) = 100 kHz 
is 50 kHz, whereas the rf frequency difference between z/ r j(5.5 mm) = 27.3 kHz and 
f r /(6 mm) = 25 kHz is only 2.3 kHz. A fine digitization is required for accuracy at small 
r that is then not necessary at large r. Method A can, therefore, be somewhat inefficient. 

Discretization method B based on equidistant steps Ai> r f : 

Writing r = r max u rf (r max ) / v rf (r) from Eq.Q gives dr = -r max v rf (r max ) / v% f du rf , and 
the integral over the radius r in Eq. ^ can be replaced by an integral over the rf amplitude 

V r f. 

u Ar \v Ar ■ \ rv r f{r min ) i 

Meff/jA t t\ max ) u r f \' mtn) I hit/ ,\ i /«\ 

n (f) = — p — _ /v / — M(z/ r/ ,t)dz/ r/ . (7) 

VrfVmin) Vrfymax) Ju rf (r max ) U rf 
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M e -^(t) can be approximated by sampling the rf amplitude v r f at N equidistant values 
is rfjk (k = 1,2,..., N) separated by Au rf = \v rf {r min ) - u rf (r max ) ] / N, with v rf}1 = 
v r f{r max ) + Az/ r/ /2 and v rf>N = i/ r /(r min ) - Av rf /2, to give 

1 - 

M e// (t)^-^ g k M(u rf , h ,t) (8) 

k=l 

with the weighting factors 

9k = -2 • (yj 

u rf,k 

In contrast to method A, method B samples the rf amplitude in equidistant steps, albeit 
with non-uniform weights for the different rf amplitudes. According to Eq.Q, the relative 
weight g k for sample k is inversely proportional to the square of the rf amplitude. Hence, 
samples with low rf amplitude contribute most to the effective magnetization vector and 
to the detected signal. For the example considered in method A, with u r f(r max ) = 25 kHz 
and u rf (r min ) = 150 kHz, the weight g{r max ) = u rf (r min ) /v rf {r max ) is equal to 6, whereas 
g{r m in) = v r f(r max )/v r f(r m i n ) is only 1/6. 



Appendix B: Optimal rectangular on-resonance exci- 
tation pulse 

For the toroid probe employed in the experiments, r min = 1 mm, r max = 6 mm, the limit 
for the maximum rf amplitude available at r max is 25 kHz, which increases to 150 kHz at 
r m in- Under these conditions, what is the optimal hard pulse length T opt that maximizes 
[M e ff (T)] x in Eq.(|3|? If we (incorrectly) choose a linear weighting of the rf amplitudes, for 
example, the duration of an ideal rectangular 90° pulse with an rf amplitude corresponding 
to \{i'rf{j' max ) + v r f(r m i n )]/2 = 87.5 kHz would be 2.86 fis. However, the non-linear 
weighting given in Eq.(|9]) emphasizes the contribution of smaller rf amplitudes, giving a 
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longer value for T^. 

A rectangular pulse of amplitude u r f, duration T, and phase y rotates the thermal equi- 
librium magnetization on resonance by an angle 27rz/ r y T about the y-axis. The opti- 
mal length of an on-resonance rectangular pulse can therefore be determined numerically 
by substituting M x (i> r f± ,t) = sin^nUrf^t) in Eq.([8]), with v r j{r max ) = 25 kHz, and 
z/ r/( r mm) = 150 kHz. We find T opt = 3.74 /is. The magnetization M x excited by this 
pulse is plotted in Fig. [8] It produces a signal that is 74% of the signal obtained by an 
ideal 90° pulse of the same length (u r f = l/(4T opt ) = 66.8 kHz) in the absence of any rf 
inhomogeneity. The performance of an optimized pulse with T = 25.62 /is, corresponding 
to the composite pulse of the same length in Fig.[TJ is plotted in Fig.|8]B for comparison. 




Figure 8: (A) The x-magnetization excited by a rectangular pulse of length T opt applied in the toroid 
probe is plotted as a function of (normalized) resonance offset v ff/v™f ax and rf scale factor r max /r. (B) 
equivalent plot for an optimized pulse at T = 25.62 /is from Fig[T] 



6 Supplementary Material 



Pulse shapes for selected pulses from the text are given as Bruker pulse shape files in text 
format. 
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3.2_CompositePPzx_4PhaseSteps 
3.2_CompositePPzx_5PhaseSteps 
3.2_CompositePPzx_6PhaseSteps 

Composite PP pulses designed according to (]22l) . Cf. Fig. 1 and 2. 
3.2_OptimizedPPzx_25point62us 

Optimized PP pulses performing the transformation M_z — > M_x. Cf. Fig. 1 and 2. 

3.3_BIR-4_50us 
3.3_OptimizedUR90_50us 

Optimized OP-BIR-4 and UR 90°_y-pulses. Cf. Fig. 5. 

3.3_OptimizedUR9(L70us 

Optimized UR 90°-pulses. Cf. Fig. 4 A. 

3.3_CompositeUR180_9PhaseSteps 
3.3_CompositeUR18CLllPhaseSteps 
3.3_CompositeUR180_13PhaseSteps 
3.3_OptimizedUR180_65us 

Composite and optimized UR 90°-pulses. Cf. Fig. 4 B. 
3.4_OptimizedUR9(L357us 

Optimized UR 90°-pulse. This pulse was used to produce the spectra shown in Fig. 7 B 
and D. 

3.5 Phase-modulated constant amplitude (v^ ax = 25 kHz) composite pulses based on the 
scheme of Ref . (|22|) : 

PP 90 x . 

Four phase steps ((f) = 225.1°, At = 2.946 ^s) 

(At)^(2A£)0_ 12 o° (2At)0_ 18O ° (2Ai)^_2io° 
Five phase steps (0 = 232°, At = 2.847 /is) 
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(At)^(2At)0_i2o° (2At)0_i 8O ° (2At) <} < ) _2io° (2At) <} i ) _225° 
Six phase steps (0 = 234.2°, At = 2.879 fjs) 

(At) ( / ) (2At)^_i20° (2At)^,_i80° (2At)^„210° (2At)^_225° (2At)<;i,_232.5 

PP 180^ 

Four phase steps (At = 1.314 fis) 

(At) o (2At)_ 120 ° (2At)_ 180 o (2At)_ 2 io° (4At)_ 225 ° (2At)_ 2 io° (2At)_ 180 o (2At)_ 120 ° (At) o 
Five phase steps (At = 1.137 /is) 

(At) o (2At)_ 120 ° (2At)_ 180 o (2At)_ 2 io° (2At)_ 225 ° (4At)_ 23 2. 5 ° (2At)_ 225 ° (2At)_ 2 io°- 
(2At)_ 180 o(2At)_ 120 o (At) ° 
Six phase steps (At = 0.966 /is) 

(At) o(2At)_ 120 o (2At)_ 180 o(2At)_2io° (2At)_ 2 25°(2At)_ 2 3 2 .5° (4At)_236.25°(2At)_ 2 32. 5 ° 
(2At)_ 225 ° (2At)_2io° (2At)_ 180 o (2At)_ 120 ° (Af )o° 
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